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i^ I An algebraic model in terms of a local harmonic boson realization was recently 
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Abstract 



proposed to study molecular vibrational spectra [Zhong-Qi Ma et al., Phys. 
Rev. A 53, 2173 (1996)]. Because of the local nature of the bosons the model 
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f^ ' has to deal with spurious degrees of freedom. An approach to eliminate the 

(N' 

^D ' latter from both the Hamiltonian and the basis was suggested. We show that 

t^ ■ 

- ' this procedure does not remove all spurious components from the Hamiltonian 

and leads to a restricted set of interactions. We then propose a scheme in 
(—] I which the physical Hamiltonian can be systematically constructed up to any 



order without the need of imposing conditions on its matrix elements. In 



^ ' addition, we show that this scheme corresponds to the harmonic limit of a 



symmetry adapted algebraic approach based on U{2) algebras. 
PACS numbers: 33.20.Tp, 03.65.Fd 



I. INTRODUCTION 



In a recent paper 0, Zhong-Qi Ma et al. proposed a model to describe the vibrational 
spectra of tetrahedral molecules and applied it to the methane molecule. The authors assign 
to each of the ten interatomic bonds a set of bosonic operators {al,ai}, i = 1, 2, . . . , 10 , 
which satisfy the standard boson commutation relations 
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The first four operators (z = 1, . . . , 4) are chosen to correspond to the C-H bonds (stretching 
vibrations), while the last six (« = 5, . . . , 10) represent the H-H bonds (bending vibrations). 
The Hamiltonian is then expanded in terms of these operators. For methane they proposed 
the following Xi invariant Hamiltonian 



H = Hi + H2s + H2b , 



(2) 



where Hi is a one-body contribution 
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Hi = es"^ altti + Cfo ^ altti + Ai ^ ajaj + A2 ^ aluj 

4 = 1 i=5 *7^i = l *7^J=5 

+A3 [(agOs + alag + a\aio) + h.c.J 
+A4 aj(a5 + ttQ + aj — ag — ag — aio) 

+02(05 — Og — 07 — Og + Og + aio) 

+a3(— as + ag — ay + 08 — Og + Oio) 



+0^( — 05 — Oe + Oy + Os + Og — Oiq) + h.C. 
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and H2X corresponds to two-body anharmonic contributions to the stretching (x 
bending {x = b) vibrations 
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s) and 



H2s = -Xs J2 N.iN, - 1) = -Xs E alaja^ai , 
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H2b = -Xb J2 NiiN, - 1) = -X, J2 <ah^a^ ■ 



(4) 
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In Eq. (0) we introduced e^ = uJs — '2,Xs and e^ = cj;, — 2X5. The number operator iVj is defined 
as Ni = altti. The Hamiltonian conserves the total number of quanta A^ = Y,l=i ^i- Since the 
methane molecule has nine vibrational degrees of freedom, while the operators {aj, aj} give 
rise to ten degrees of freedom, spurious interactions may occur in H. Their identification, 
however, is not obvious when the Hamiltonian is written in the form of Eq. (0)-(|^). In ||l| the 
authors proceed by identifying the spurious states as those containing the fully symmetric, 
one-quantum bending state i^^Ai, 100000) as a factor. With a combination of projection 
and coupling techniques and an orthogonalization procedure, the physical space up to three 
quanta is generated. Next it is shown that the physical Hamiltonian does not couple the 
spurious states with the physical states as long as the conditions 

A5 = , Xb = , (5) 

are imposed on the Hamiltonian of Eq. (^. Although these conditions indeed lead to 
the correct results, they are too restrictive and, furthermore, do not remove all spurious 
interactions. 

II. SPURIOUS INTERACTIONS 

In this section, we show that the condition of Eq. (^ is not sufficient to remove all 
spurious interactions from the Hamiltonian. We define a new set of harmonic bosons with 
well defined tensorial properties under the tetrahedral Td group by projecting the local boson 

operators. For the stretching modes we have 

4 

with r = Ai, F2 and for the bending modes 

10 
&r,7 = E A"'" a\ , (7) 

with r = Ai, i?, F2, and similar expressions for the annihilation operators. The coefficients 
Oj '''' and /3j ''*^ are the same as those appearing in the projected wave functions of Eq. (3.1-2) 



in [jT|. In the framework of this alternative scheme, both the Hamiltonian and the basis 
states are constructed by repeated couphng of these normal bosons. In this way, we obtain 
automatically an orthonormal symmetry adapted basis with a well defined number of quanta 
in each one of the fundamental vibrational modes. For example, with the usual notation 
for the normal modes of methane (z/i, 1^2, u^, u^), we find for two quanta the following basis 
states with symmetry Ai 

2,0,0,0) = -^(44)^^10), 
0,2,0,0) = 4^(fel,6l,)^^|0), 
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0,0,2,0) = -^(44)^^ 

0,0,0,2) = i=(44)^^|0) 
0,0,1,1) = (46l.J^^|0), 



(8) 



Here (. . .)^ indicates coupling of two quanta to symmetry F. These basis states are orthog- 
onal, since they correspond to eigenfunctions associated with different eigenvalues of the 
number operators 



^r, = J2 5r,7^r,7 , ^r, = ^ &r,^&r, 
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(9) 



This symmetry adapted procedure provides a method to remove the spurious interactions 
from the Hamiltonian as well as the spurious states from the basis. It consists in deleting the 
tensors which are associated with the spurious bending mode, b\_^ and 6^^ , in the construction 
of physical operators and basis states. This allows an exact elimination of the spurious states 
without the need to impose additional conditions on matrix elements. 

The Hamiltonian of Eq. (^ does include spurious interactions. This can be seen by 
rewriting H in terms of the normal bosons of Eqs. (^ and (|^). For the one-body Hamiltonian 
of Eq. (^) we find 

Hi = {es + 3Ai) fiA^^ + (e^ - Ai) fip:,, 

+ {eb + 5X2 + A3) uau + (efe - A2 + A3) fiE, + (efe - A2 - A3) hp^^ 



+2v^A4 lisibpy'+h.c] +2V6X 



'>2"-F2 



[46^J^i+h.c. 



(10) 



There are two terms that contain the spurious Ai bending mode. This shows that, in 
order to remove the spurious interactions from Hi, the condition A5 = of Eq. (^) is not 
sufficient. In the present symmetry adapted analysis we find two conditions: A5 = and 
efe + 5A2 + A3 = . We remark that, since the expectation value of tiaj^^ in the physical basis 
vanishes, the omission of the second condition has no further consequences. 

A similar analysis can be carried out for the interactions in H2S and H2b- The opera- 
tors in Eq. @ represent strong local interactions, and their representation in terms of the 
normal operators of Eqs. (|^) and (|^ acquire a complicated form which involves operators 
that transfer quanta between different modes. The Hamiltonian H2S only depends on the 
stretching modes and hence does not contain spurious interactions. However, the bending 
contribution i/26 contains terms that transfer quanta between spurious and physical states. 
This can be seen by calculating its matrix elements in the basis states with Ai symmetry 
that have two quanta in the bending vibrations 



|r) = ^(«)^^|o), 



(11) 



with r = Ai, E and F2, respectively. The states with T = E, F2 belong to the physical basis 
of Eq. (§), whereas the state with F = y4i is spurious. In this basis the matrix elements of 
H2b are 
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(12) 



The nonvanishing matrix elements in the first column and the first row show that H2b con- 
tains terms which involve the spurious Ai bending mode. For this reason the parameter 
Xb in Eq. (^) must vanish in order to remove the spurious contributions from the Hamilto- 
nian, in agreement with the result obtained in |I[]. However, this condition also removes the 
physical part of H2b that couples the states with T = E and F2. 



In summary, in this section we have shown that in order to remove the spurious interac- 
tions from the Hamiltonian of |jl| three conditions have to be fulfilled 

As = , Xb = , e^ + 5A2 + A3 = . (13) 

In the calculation of the vibrational spectrum of methane presented in [^, no stretching- 
bending interactions were included {Xb = 0). The spurious components in the Hamiltonian 
were removed by requiring that its expectation value in the spurious state with one quantum 
in the Ai bending mode vanishes. This method only leads to exact results in the harmonic 
limit in which it reduces to the last condition of Eq. (|T3p. 

An additional comment concerns the basis used in ^. The basis functions in general 
do not carry the quantum numbers that correspond to the number of quanta in each of the 
normal vibrational modes. This represents a serious disadvantage in the process of identify- 
ing the calculated energies with the experimental ones, which traditionally are reported in 
terms of these labels. This problem was overcome in p[ by constucting the physical basis 
by means of the %i Clebsch-Gordan coefficients. We note, however, that in Eq. (3.12) of 
the normalization factor 



|n)|m) = Y[ 



( I Ml 1/2 



rif^lmi^l 



n + m) , (14) 



was missing. When this factor is included, the coupling procedure provides an orthonormal 
symmetry-adapted basis for the physical space with the desired labels (z/i, 1^2,^3, ^i)- 

III. SYMMETRY ADAPTED APPROACH 

At this point we would like to mention that recently the ^(2) model used in 0] has 

been further developed |§|-^ by elaborating a tensorial formalism analogous to the one 

suggested above. In this symmetry adapted approach the local f/j(2) algebras {Gj} = 

{Ni, J^^i, J_^i, Jo,j} are projected to % tensors 
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T',:, = E «r'" J,. , Ti^ = E A"'" I. > (15) 



where /i = ±,0. The coefficients have the same meaning as in Eqs. (j^) and (J^). The 
Tct invariant Hamiltonian is constructed by repeated couphngs of these tensors to a total 
symmetry Ai. For example, the interactions that are at most quadratic in the generators 
and conserve the total number of quanta are given by 
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V.-^EtS^T-V (16) 

Here T = Ai, F2 for stretching vibrations (a; = s) and T = E, F2 for bending vibrations 
{x = h). The harmonic limit of these interactions is found by taking the limit Ni -^ 00, so 
that [|,| 

lim J!- = tti , lini jl. = a] , lini = 1 , (17) 



where Oj and a] are the boson operators used in [1[]. In this harmonic limit the interactions 
of Eq. ( ]T6D take the form 
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lim Tir^ = 77 H (sr,7 ^r,7 + sr,7 « 
lim Tir, = - 5^ (^r.^, &r,7 + f'r,7 ^r 



lim Vr. = , 
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n}}%Ioo ^'^=2^ (^k,7 ^^2,7 + SF2,^ b^F2,yJ > 



lim Vsb = ■ (18) 

which corresponds to the interactions appearing in Eq. (|10D, with the exception of the spu- 
rious interaction which, by construction, does not arise in this approach 0. The symmetry 
adapted basis states are constructed as in 0, but now taking into account the appropriate 



normalization coefficients of Eq. (0). The application of this model to the vibrational spec- 
trum of methane yields a fit which is an order of magnitude more accurate than those of |]1| 
and §. 

IV. SUMMARY AND CONCLUSIONS 

Summarizing, we have analyzed the boson realization model presented by Zhong-Qi Ma 
et al. in |jl| and found the following results: 

(i) The approach presented in [Q does not remove all spurious interactions from the 
Hamiltonian. Moreover, the condition Xb = eliminates, in addition to spurious interac- 
tions, also some physical interaction terms. 

(ii) Although the method proposed to remove the spurious states gives the correct results, 
it becomes increasingly difficult to apply when additional higher order interactions have to 
be included in the Hamiltonian. 

(iii) In general, the basis in which the Hamiltonian is diagonalized does not carry the 
quantum numbers associated to the number of quanta in each fundamental vibrational 
mode. This problem can be solved either by using the T^ Clebsch-Gordan coefficients or by 
diagonalizing the number operators of Eq. @. 

(iv) The formulation in terms of tensor operators provides a natural way to eliminate 
the spurious interactions from the Hamiltonian. Besides, we have shown that the approach 
of [jll corresponds to the harmonic limit of a more general method based on coupled U{2) 
algebras P-§]. 
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